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The metric recommanded by the IAU2000 resolutions allows propagation of light calculations at 
the c -3 level in the general relativity framework. In a recent paper [I], motivated by forthcoming 
space experiments involving propagation of light in the Solar System (ASTROD, GAIA, LATOR, 
ODYSSEY, SAGAS, SIM, TIPO, ...), we have proposed an extention of the IAU metric equations 
at the c -4 level. This has been made in the general relativity framework. However, scalar-tensor 
theories may induce corrections numerically comparable to the c -4 general relativistic terms. Ac- 
cordingly, one proposes in this paper an extension of [I] to the scalar-tensor case. The case of a 
strongly hierarchized system (such as the Solar system) is emphasized. In this case, an explicit 
metric solution is proposed. 
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I. INTRODUCTION 

Forthcoming space missions and missions in project 
(like ASTROD, TIPO or LATOR, see [T]) will require 
distance measurements at millimetric level in the Solar 
System. This corresponds to time transfer at the preci- 
sion 10 _u s. As argued in p], this requires a complete 
Solar System metric at the c -4 level, in order to describe 
the laser links involved in such experiments. This has 
been proposed in the framework of the General Relativity 
(GR) theory in [1], leading to an appropriate extention 
of the metric recommanded by the IAU2000 resolution 

The relative amplitude of the relativistic effects is of 
the order of e = GM/rc 2 , where M and r are some char- 
acteristic mass and distance. The so-called first order 
terms are of order e, the second order terms of order e 2 . 
In the inner Solar System, e is typically of the order 10~ 8 , 
and can be sensitively greater (10~ 7 ) for photons entering 
well inside Mercury's orbit, and can even be as large as 
10~ 6 for photons grazzing the Sun. On the other hand, 
GR faces a lot of difficulties both at the solar system level 
(Pioneer anomaly, fly-by anomaly) and at the cosmolog- 
ical level (universe expansion accelarate). Hence a new 
surge of interest in alternative gravity theories. Among 
these, the scalar-tensor (ST) theories deserve a particu- 
lar interest, since the gravitational sector of fundamental 
theories, like string or Kaluza-Klein theories, turns out 
to be described by a metric tensor plus a scalar field [3J H] 
(Brans-Dicke or not [5]). 

Alternative gravity metrics diverge from the GR one at 
the c -2 level. This divergence is quantified by the Post- 
Newtonian (PN) 7 factor (— 1 in GR) entering the c -2 
term in the space-space components of the metric tensor. 
From the present observations, |7 — 1| can at best reach 



values of the order of 10 -5 [7]. Besides, some theoretical 
considerations strongly suggest 7 could be driven from 
any "initial" value to unity by the cosmological expan- 
sion (more precisely, ST theories are driven to GR, as 
soon as these theories fulfil some (not very constraining) 
conditions), and also that the I — 7 current value could 
be of order 10" 7 or ICT 8 [HIS]. On the other hand, c~ 4 
terms are typically 10 6 to 10 8 smaller than c~ 2 terms. 
Hence, as much as c~ 4 space-space metric terms have to 
be taken into account in light propagation problems, it 
is necessary to include also a possible divergence from 
GR in the proposed metric, at least in the c -2 terms, for 
numerical coherence. Hence it is required to upgrade the 
(extended version [T] of) IAU2000 metric to the ST case. 

In principle, this would call for new definitions of mul- 
tipolar moments at the c -2 level. But other publications 
porposcd c~ 2 mutlipolar moments in ST theories |10j or 
in a parametrized post-newtonian framework [11] . Thus 
we do not discuss this point in this paper and focus only 
on the c metric side of the problem. 

In section |Hj one defines a terminology relevant to the 
considered problem. Since the Einstein conformal rep- 
resentation plays a central role in the approach followed 



in this paper, the section [Hi] is devoted to the confor- 
mal link between the representations of the ST theories, 
and the related notations we will use. Section HVl is ded- 
icated to the derivation of the ST field equations up to 



O (c 5 ) terms. In section JvJ these field equations are 
rewritten up to 0{c~ 5 ,uJq c -4 ) terms, for applications 
taking explicitely our present knowledge on ST theories 
into account. Finally, considering applications to So- 
lar System-like systems, one defines strongly hierarchized 
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system in section VI This is made by defining a quantity 
fi that quantifies how much the system is gravitationally 
dominated by its most massive body. Accordingly, the 
field equations are rewritten up to O (c~ 5 , Wq^c -4 , /ic~ 4 ) 
terms. In this case, the explicit solution that can be used 
in relevant applications is written in well-suited coordi- 
nates. 
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II. TERMINOLOGY : DEFINITION OF THE 
PN/BM AND PN/RM METRICS 

The PN approximation is based on the assumption of 
a weak gravitational field and weak velocities for both 
the sources and the (test) body (ie. velocities of the or- 
der y/GM] r or less, M being some caracteristic mass of 
the system). It formally consists in looking for solutions 
under the form of an expansion in powers of 1/c. The 
usually so-called nPN order terms in the metric, lead- 
ing to c~ 2n terms in the equation of motion of a body 
describing a bounded orbit, are terms of orders c~ 2n ~ 2 
in goo, c~ 2n ~ 1 in g^ and c -2 ™ in gij. In this case, the 
Ricci tensor components have to be developped the same 
way as (18) in pQ. In this paper, a metric developped 
this way will be refered as the nPN/BM metric (BM 
meaning "Bounded Motion" for test particles). It is 
particularly well-adapted for studying bounded motions 
in systems made by non-relativistic massive bodies, as 
the Solar System is. 

However, since we are interested in the propagation of 
light, we are lead to relax the hypothesis on the velocity 
of the test particle whose motion is considered. Of course, 
this doesn't change the metric, but the terms to be con- 
sidered in the metric components are not the same as in 
the PN/BM problem. Indeed, the terms leading to c~ 2 " 
terms in the equation of motion of a test particle moving 
with relativistic velocity are terms of order c -2 ™ in all 
the components g a p (ie. in both goo, 9oi and g^). In this 
case, the Ricci tensor components have to be developped 
the same way as (19) in [T]. In this paper, a metric de- 
velopped this way will be refered as the nPN/RM metric 
(RM meaning "Relativistic Motion" for test particles). 
A PN/RM metric is particularly well-suited for study- 
ing relativistic motions of test bodies (for instance, the 
propagation of light) in systems made by non-relativistic 
bodies, as the Solar System is. 

PN/BM and PN/RM orders are illustrated in figure [l] 

The present paper deals with the PN/RM problem 
since we are concerned in propagation of light. 



III. THE ST THEORIES IN EINSTEIN VS 
JORDAN REPRESENTATIONS 

The Jordan representation of the ST theory is de- 
scribed by the action 



S = 




16ttG 

+ / d i Xyf^gL NG {^,g IJ , v ) 



<]>/,> _ ^Ig^d^dp^ 



<I> 



(1) 



In this representation, the gravitational sector of the the- 
ory is described by the Jordan metric g a p and the scalar 
field while the non-gravitational fields are symbol- 
ically represented by ^ . The scalar field couples di- 
rectly with the metric, leading to rather complicated 
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FIG. 1: General scheme of orders taken into account in 
PN/BM and PN/RM metrics 



field equations. Besides, the kinetic term associated to 
the scalar field doesn't have the standard form, involv- 
ing a scalar field function w, the function characterizing 
the ST theory we are dealing with. On the other hand, 
the non-gravitational lagrangian Lng doesn't depend on 
the scalar field, leading to simple equations of motion 
(V a T al3 — 0), with the nice consequence that the weak 
equivalence principle applies in this representation of the 
theory. 

To overcome the just mentionned drawbacks, one could 
be tempted to resort to a dependent variables change 

(g a p,$) — > {9 a f},v) 
chosen in such a way that the action transforms into 



S 



16itG 



d xy— g 



R - 2g a P d a ipdpip 



+ J d i x^/^g~L NG 



(2) 



Lng depending on g and tp in a way to be precised 

later (i? and g correspond to R and g, but with g a p re- 
placed by g^a). Since the scalar field doesn't couple with 
the metric, (pi) is refered to as the Einstein representation 
of the theory. 

The form ([2]) is achieved by considering a conformal 
transformation of the metric 



g af} =A(ip) 2 g a0 . 



From the induced transformation of the Ricci scalar 



(3) 
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and up to a divergence term, ^ turns into 



# A 2 \Ja) j9 aP daAd p A] 

d A X\/—gA~ i LNG ■ 



(4) 



Comparing Q with ([I]) suggests : 

- the link between the Jordan and Einstein represen- 
tations of the scalar field 



(5) 



the identification 



2 f dtp 
A?(dA 



g ap d a Ad p A; 

. (6) 

- the Einstein representation of the non gravitational 
lagrangian 

Lng = ^M 4 ire(*,v) 

= A (<p) 4 L NG A (cp) 2 g^ . 

§5§ and ^ lead to the link between the functions u (<&) 
and A (ip) (equivalently characterizing the considered ST 
theory) 



(2lu + 3) 



dln\A\ 

dip 



(7) 



([7]) requires u) > —3/2. This results from the fact one has 
imposed the sign of the scalar kinetic energy term in Ein- 
stein representation ^ in order to ensure the dynamical 
stability of the theory [T3"] . 

The link between the stress tensor components in the 
two representations follows from the general stress tensor 
definition 



a/3 
a/3 



S I d xy/^gL NG : 
S I d i x^J^L NG ee - 



a/3 



d x^gT ap _ 
d A xy/^ft a pSg a 



Since \/ —gL^a — \/~gLNGy and since the scalar field is 
not varied in this metric variation process (no ambiguity 
since the two versions <p and $ of the scalar field are 
related in the non metric dependent way (|5|), it directly 
turns that 

T a p = A 2 T a p. 
For the mixed and contravariant components, it follows 



T a = A^T^ ( =>> f = A i T), f 



a/3 



j^rpa/3 



the indexes being raised/lowered by the metric involved 
in the corresponding representation. 

Eliminating A 

It is clear A can be eliminated between the two repre- 
sentations of the scalar field using ^. This way, the con- 
sidered ST theory is represented by the function $ (ip) in 
Einstein representation. From ^ , this function is linked 
to lu ($) by 



(2oj + 3) 



dln<f> 

dip 



(8) 



The conformal transformation ([3|, the link between the 
non gravitational lagrangian and stress tensor represen- 
tations now write 



g a /3 = ®9a/3 



L NG = <P 2 Ljv G (*,$ l g ) 



T a /3 = $ 1 T a p 
-J 



(9) 



T 



a/3 



$-3 T a/3 



(10) 



As it turns from equation ([8]), <p is defined up to a sign 
and an additive constant. When needed in the following, 
the sign will be fixed by the choice 



-eHn$ 

dip 



(11) 



IV. THE FIELD EQUATIONS 

A. Using the Einstein representation gravitational 
held variables 

From pj), the field equations can be written 

R = ^(f^-ifp^+MV^ 



On 



'a (y~^d ip) 



2-kG - r^dhx® 
-^ TV ~ 9 ^p- 



(12) 



the function $ (ip) characterizing the ST theory explic- 
itly entering the scalar field equation. As usual in PN 
approximation, let us write the scalar field as 



(2) (4) 



(13) 



(2) (4) 

where ipo is constant and ip and p arc zcroth order 
terms. (Remark that, since ip is defined up to an addi- 
tive constant, it is not restrictive to set ipo — 0.) As a 
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consequence, it turns out that, under the standard PN 

- ij . 

assumptions, R = O (c J , so that the Strong Spatial 
Isotropy Condition (SSIC) [H] applies in this represen- 
tation. It is then possible to choose a coordinate system 
in which the Einstein metric takes the following form at 
the 2PN/RM level [1] 



ffoo = 
9m = 

9 in = 



2w 



2w 2 



4wi 

^,3 



1 



0(c 
2w 

72" 



2w 2 



4t, 
c 4 



Oic 



The scalar field function develops as 

* ft,) = i + + i (*$ + l -Kf) + o ( c - 5 ) 

(15) 

where $ (= 1), & and $ ' stand for the values of $ 
and its derivatives at <po. Putting $o = 1 is not restric- 
tive since <& enters ( 12 ) throught its logarithm derivative. 
Setting 



1 



(T* 



oo 



(16) 



cr 



( 



Jfefc 



(which, from standard PN asumptions, are c° order quan- 
tities) and using (10 1, the (00), (Oz) and (ij) field equa- 
tions (Il2| lead respectively to 



Aw- 



1 



(3d tt w + AdtkWk) , 
cr c 



3 $; ( ^Aic = -47rGcr+0 (c~ 3 ) 



Awi — dikWk — daw — —AirGa 1 + O (c 



(17) 
(18) 



Remark that, in contrast to the GR case, the (00) equa- 



tion is not linear, because of the c <p Aw term. The 
scalar field equation also contains a non-linear c -2 term. 



Now, combining (171 and (21) leads to 



A V 



(, 2 ) 1*/ 



(14) Accordingly, let us put 



(2) 
9 



(4) 

Hence, defining \ = f /®'o 



2c 2 

The metric field variables 



(22) 



(23) 



W, Wi 



and Ti, are now decou- 



pled from the scalar (remaining) one x- T ne system con- 
straining w, Wi and r,j now writes 

Aw+\ (3d tt w + Ad tk w k )+^<P Q wAw = -AnGa+O (c~ 3 ) 

2c z 



Awi — dikWk — duw = —AirGa 1 + O (c 



©ij (m) = -duWj - dtjWi 
1 



(24) 
(25) 

(26) 



1 



- 28 l3 (duw + d tk w k ) + AirGa 1 ' + O (cr 1 ) 
X being obtained in a second step, by solving 

A X - 2(d tt w + d t kw k ) (27) 
+ 1 ($^ 2 - $£) wAtu = -2 7 rGa' £fc + O (c^ 1 ) . 

Remark that Aw may be replaced by — AixGa in the non 
linear terms of equations (24) and (27). 



©ij (m) = diwdjW - d t (diWj + d 3 Wi) 

- 2Sijd t (d t w + d k w k ) (19) 



+4nGa ij + d-Xpd-Xp + O (c" 1 ) 



where O^- is defined, as in pQ, by 

©ij (tm) = dikTjk + djknk — At. 13 — d l3 Tkk- 
The scalar field equation gives 



A (2) 1 f S. (2) , A (4) 

A ip + [ -d tt ip + A ip 



1 ( . $ '\ ( 2 ) A ( 2 ) 



= 2 rC ,'<!>;, ( a + — | +0 (,- 



-3 s ! 



(20) 



(21) 



B. Back to Jordan representation 

Let us use the function uj ($) and its derivative u/ (<£>) 
instead of $' (y>) and (y>). One finds, using ([TT]) 



V2wo + 3 



2w + 3 



1 - 



2w + 3 



4^o 



" ~ (2 Wo + 3f 



One now goes back to Jordan representation using (pj), 
with, from (15) and (23), 



2w 



c 2 (2u + 3) 



(28) 



(2w + 3) 



2wn 



2wn 



w —Ax 



+ 0(c" 5 ). 



5 



Now let us put 



and, for the matter part of the equations 



7 

G eff 

and let us define 



w + 2 

1 + < 

(2uj + 3) (2uj q + 4) 2 

2cj + 4 



£ = 



1 



2wn 



-G 



(U,Ui,Uij,P) = ^3 { w i w u T ij,X) 
and the related quantities (W, Wi, Wij) by 

W = C/ + (l- 7 )4 
<r 

Wi = Ui 

= Uij-il-iiPSa. 

Using (14 1 and (281, one gets the Jordan metric 



5oo 
9oi 



1 2W 2W 2 „ i -5^ 
-( 7 + l)^ + 0(c- 5 ) 



(29) 



r l\ 2^ /j fl ,n 2W 2 

<fc = 1 + 7"^ + (7 2 +/? - 1) —r 



(7+1)^+0^ 



where the functions (W, Wi, Wij,P) satisfy the following 
field equations - after some algebra from ( 24|27 l and ( 16 ) 



DW + 1 + 2/ ! 3l WAW + % (1 + 7 ) d t J 
c c 

= -^G eff Z + O (c- 3 ) 

AW t - d t J = -4 7 rG e //S i + O (cT 2 ) 

AWij + d t WdjW + 2(1-/3) S^WAW 

- diJj - djJi - 2-ySijdtJ = -47rG e// £« + O (c 

AP + 2^^WAW ~ 2d t J 



1-7 

= _47rGe ^3 7 ^I + ° (r 



(30) 



In ( 30 1 , one has set 



J = d t U + d k U k (31) 

= 9 t W + d k W k + O (c~ 2 ) 

Ji = d k U lk - -diUkk + d t Ui 

= d k W ik - \diW kk + d t Wi - ] ^r 1 d t P (32) 



(T°° + 7 T fcfe ) 



E 1 = -T u 
c 



S« = T» J - tT**^ ( 



(3 7 - l)T fcA: ) 



Let us remark that the quantity 2 (/3 — 1) / (1 — 7 ) (= 
w (2wo + 3) _ (2w +4)~ ) is not diverging, even if 7 
is (arbitrarily) close to unity. Besides, no new PN pa- 
rameter appears neither in the c~ 4 space-space part of 
the metric nor in the corresponding field equations, as 
stressed in [5]. 

This form is relevant in all sufficiently weak gravita- 
tional field, even in systems where the ST theory is not 
very close to GR, i.e. the PN parameters 7 and j3 are not 
necessarily very close to unity. A priori, this may occur 
even if 7 and j3 are close to unity in some (other) regions 
of the universe, as in the Solar System, as soon as the ST 
theory is not (in some sense) close to the Brans-Dicke one 
(in Brans-Dicke gravity, U) doesn't depend on the scalar 
field, so that it as the same value in all the space-time 
regions of the universe). 



C. Harmonic gauges 

Since the use of the harmonic gauge (HG) is recom- 
mended by the IAU, let us consider the field equations 
in this gauge. Of course one has to specify the represen- 
tation in which the HG is prescribed. The Jordan HG 
condition reads 

9 1 afl — U 

and it leads to, for the space (a = A;) component 

( 7 -l)d k U = 0(c- 2 ). (33) 

As expected from known results in GR [14j . this condi- 
tion reduces to a triviality in the case 7 = 1. On the other 
hand, ( 33 ) with 7^ 1 shows that the coordinate system 



in which the metric takes the (Jordan) form ( 29 ) , corre- 
sponding to SSIC in Einstein representation, doesn't en- 
compass (Jordan) harmonic coordinates in the ST case. 
In other terms, (Jordan) harmonic coordinates are in- 
compatible with the SSIC in Einstein representation. 

One could rather choose to impose the HG condition 
on the metric in Einstein representation 

9 r a /3 = 



since the Einstein metric ( |14[ ) satisfies the SSIC. That 
means one imposes to, tOj and Tij to satisfy, in addition 
to the field equations, the relations 

d t w + d k w k = o(c- 2 ) 
d k r ik - -diT kk + d t Wi = O (c -1 ) . 



G 



Translated in terms of (U, Ui, Uij), this takes exactly the 
same form, i.e., using ( 3l]|32 ) 

J = 0(c' 2 ) , J^O^ 1 ). (34) 

It turns out this corresponds to the Nutku gauge con- 
straints ([10], [H]), meaning that imposing the HG in the 
Einstein representation is equivalent to impose the Nutku 
gauge in the Jordan representation. Using (34), the three 



first equations of (30) take the reduced form 



nw + 1 + 2 ?-^wAw 

cr 

= -47rG e// I] + 0(c- 3 ) 
AW, = -47rG e// S 4 + O (c~ 2 ) 
AW i:j + d.WdjW + 2(1-/3) SyWAW 

= -inG eff J: lj +0(c- 1 ) 



(35) 
(36) 

(37) 



while the equation corresponding to the scalar degree of 
freedom and the harmonic constraints read 

5 — 1 Y, kk 
AP + 2- WAW — —AivGeff- + 0(c~$8) 



1-7 " w J 37 - 1 

d t w + d k w k = o(c- 2 ) 
d k w ik - + dm - x -^- 



diP = O i 



(39) 
- 1 )(40) 



or equivalently (after elimination of the scalar field P) 
d t W + d k W k = O (c- 2 ) (41) 
d ik W ik - \AW kk + d ti Wi + 03-1) WAW 

= -2vG eff ^^ kk +0(c- 1 ) (42) 

d ik W jk + d ti Wj = d jk W ik + d tj Wi + O (c" 1 ) .(43) 

The last equation refers to the fact that diP (given by 
(40)) is a gradient. 



One could note these equations are coherent with 
1.5PN/BM equations assumed in [TT] . 



V. SIMPLIFIED FORMS THAT CAN BE USED 
CONSIDERING PRESENT CONSTRAINTS ON 
GRAVITY 



useful metric resulting from (29) reads (if uj'q is not "un- 
reasonably large") 



2W 2W 2 

4~ 



9oo 


= -1 


+ —2 


9o% 




9ij 


= 


2W 

1+7— 
c 1 



0(c- 5 



m. 



l„-4> 



where W, Wi and satisfy, from (30) 



DW + -^d t J = -AirGeffcr + O (c^ 3 , c^c" 2 ' 

AW t - diJ = -4nG eff (T l + O (c~ 2 ) 
AWij + diWdjW - d.Jj - djJi - 28^ d t J 
= -4irG eff * l i +0(c-\luv 1 ) 

and where J,- reduces to 



Ji = d k W lk - -dtW kk + d t Wi + O (uio 1 ) 



(44) 



One remarks the field equations take exactly the same 
form as the GR case p] (with G replaced by G e //). The 
only remaining reference to the scalar field is reduced to 
the presence of the 7 PN coefficient in the metric tensor. 
Related to this, the field equation on P is dropped out. 



A. Harmonic gauge 

The corresponding harmonic equations to be used 
when considering known constraints on loq reads 

DW = -AnG eff a + O (c -3 , w^V 2 ) 

AW, = -4irG ef f a l + O (cr 2 ) (45) 

AWij + diWdjW = -4TrG eff a t: > + O (c _1 , uJq 1 ) 

with gauge conditions 



In the inner solar system, the gravitational field is such 
that 



2U 
72" 



10" 6 to 10"" 



On the other hand, from experimental/observational con- 
straints [7] 

|7- 1| ^ !0~ 5 ie. w > 10 5 . 

This means 7 — 1 (or Wq 1 ) could be considered numeri- 
cally as a c _1 (at best) order quantity. Hence, it is con- 
venient to present the metric under the form of a gener- 
alized development in both powers of c _1 and Wq" 1 . The 



d t w + d k w k = o ( c - 2 ) 



(46) 



d k W lk - -d,W kk + d t W, = O {c-\^ x ) . (47) 



VI. EXPLICIT FORM TO BE USED IN THE 
CASE OF THE SOLAR SYSTEM 

A. Strongly hierarchized systems 

Let us consider the case where the system is composed 
by bodies of masses Ma- Let us consider one of these 
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bodies, named 5, of mass Ms- Let us define the param- 
eter 



r B £ Ma - 



M s 



A^S 



Let us also point out that, since at this level the metric 
depends on time through the positions of the star and 
the planets only, the operator d t is of order 0{y). Hence, 
equations ( 45 1-( 47 ) lead to 



One defines a strongly hierarchized system as a system 
in which it is possible to choose the body S in such a way 
that 

H < 1. 

In such a system, the body S will be hereafter referred as 
" the star" , while the other bodies will be referred as the 
"planets". 

In the general relativistic A^-body problem, multipolar 
moments of a body A are defined in the coordinate sys- 
tem in which this body is, in some sense, at rest. These 
moments are affected by coordinate transforms through 
a "Lorentz-like length contraction effect". These effects 
being of order (u/c) , where u is the relative velocity be- 
tween the two frames, the induced effects in the metric 
components are of order c~ 4 . 

In strongly hierarchized systems, the velocity of the 
body S is of the order of 



V S ~ fiVB ~ Ml 



' GMs 
rs-B 



where B is the most massive planet (and rs-B the dis- 
tance between B and the star) . All the Lorentz-like con- 
traction terms have the form 



GM A v\ 



If A is a planet {A ^ S), this term is of order O (/xc" 
If A is the star (A 
ie. O 



S), this term is of order O 



VaC 



[i 2 c 4 ) . Hence all these terms are, at best, of order 



AW = -47rG e// cr + (c^ 3 ,^ 1 ^ 2 ,^ 2 ) 

AWi = -ATtGeffd 1 + O (cr 2 ,^- 1 ) (48) 

AWij + diWdjW = -AirG ef fa ij + O (c -1 ,^ 1 ^) 



with gauge conditions 



d t W + d k W k = 0(c- 2 ,fic- 1 ) (49) 
dkWi k -la i W kk = 0(c- 1 ,Uo 1 ,ii). (50) 



Related to this, defined in [5] leads to numerically 
negligible terms (see (11.4.8) in [10] for the ST version). 
B. Application to the solar system 

In the Solar system, the most massive body S is the 
Sun and one has 



Thus, it is legitimate to consider the Solar system as a 
strongly hierarchized system. Note that, at best, only 
J2, J4 and Jq planetary terms (giant planets) could have 
a significant impact on laser ranging experiments at the 
required accuracy (see [2iTf24] for giant planets' multi- 
pole moments values). Hence, taking advantage that the 
Solar multipolar terms are very weak, the solution of the 
field equations (48) - suitable for millimetric accuracy 
in propagation of light - with the harmonic constraints 
( 49p0 ) given in barycentric coordinates turns to be 



where 



2 2W 2 
.900 = -l+^[W (t,x) + W L (t,Z)}- —^ + 0(c- 5 ,uJo 1 c-\Lic-\j$c- i ) (51) 



90i = -2^ r W l (t,x)+0(c- 5 ,nc- i ,Jic- i ) (52) 
9ij = (l + J f ) + W & f )l + ) S * + 4 ^T + (c" V ~ ^ A^ 4 , J|c- 4 ) (53) 



W (t,x) = J2 W A,o> with W A ,o = G eff M ^ (54) 

3 H 2n ( k A ■ f A \ 

W L (t,x) = ^^, L , With W A ,L = -Geff^MAJt ^TPT P 2n ■ (55) 



W % {t,g) =^2wX(t,£), with W A (t,x) 

A 

W s (x) 
W l3 {x) 



G 



eff 



[r A x Sa^J 



2r\ 



M A v\ 
rA 



ff - 



Ms 
rs 



3 j^2n 
1 + X! J 2n-£^ P 2n 



1 (r M s y {x i -x i s ){xi-x> s ) 
4 [ Geff ^ 1 1 " ' 



rA 



(56) 

(57) 
(58) 



r 



where one has put 

f A (t,x) =x-x A (t) and r A (t,x) = \r A (t,x)\. 

Ma, r a , va and Sa being repectively the mass, the po- 
sition and the velocity in barycentric coordinates, and 
the total angular momentum of the body A. R A and 
j£ are the radius and the mass multipole coefficients of 
the body A. P n are the Legendre polynomials and k A 
denotes the unit vector along the local Z A axis of each 
body A. The differences with the IAU2000 metric [2J lie 
in the presence of both the PN parameter 7 and the c -4 
space-space metric term. The multipolar term Wl in <?jj 
that has been neglected in the IAU2000 metric - thanks 
to numerical considerations in the 1PN/BM case - has 
to be considered here as well. 

In a same way, the multipolar terms in the time-space 
component of the metric (goi) could also lead to measur- 
able effects. Thus, one may have to consider them. 

While the rotational term in the time-space compo- 
nent of the metric is given as the usual Lense-Thirring 
term, slight modifications (spin multipoles) can in princi- 
ple appear due to the differential rotation of the bodies. 
However, Solar seismology suggests [25] that the Sun's 
tachoclyne is at about 0.7 Sun radius. Then, the mass 
concerned by the differential rotation is of order of a few 
percent of the total mass and thus, it might not lead to 
measurable effects. But, incidentally, time transfer and 
laser ranging experiments could suggest a way to test 
our knowledges on the solar interior dynamics, indepen- 
dently of results coming from Solar seismology. Another 
point coming from the solar seismology is that the solar 
core (r < 0.2 Sun radius) may rotate faster than the ex- 
ternal layers [27H2D] . Since, the core represents a great 
amount of the total mass, this could affect the propaga- 
tion of light at a level depending on the total angular 
momentum value (Ss)- Depending on the solar inter- 
nal structure model, this may happen at the millimetric 
level. Using a simplified model, we show in appendix [A| 
how a non-rigid rotation of the Sun could also affect the 
metric. 

Let us note that Wl can be neglected for inner so- 
lar system millimetric laser ranging experiments, such as 
Mars laser ranging for instance. 

We emphasize again that neiher the ft parameter nor 
the e parameter are required (e - corresponding some- 
times to A [3DJ or 5 [ID] - being some PN parameter often 



considered in the c~ 4 space-space metric term [31 ] 132 ) ). 
This is because both the former and the latter give too 
small deviations from GR to be considered in c~ 4 So- 
lar system photon's trajectories calculations. This fact is 
known for the former from [7] and is then obvious for the 
latter since it is a function of 7 and ft in (non-massive- 
)ST theories considered here (as expressed in equation 



Finally, let us point that, because of the orders of mag- 
nitude in the Solar system problem, while we must know 
the global (ie. BCRS) 2PN/RM metric for experiments 
of interest - which is given in (|51[)-(58| - it is not neces- 
sary to have the local (planet) 2PN/RM metrics, insofar 
such metrics are needed. This means that the usual lo- 
cal 1.5PN/BM metrics ([H] for GR and [TO] for ST) are 
accurate enough for overlaping local coordinates charts 
with the global one. Hence, it is not needed to develop a 
complete theory of reference frames at the 2PN/RM level 
as long as one is interrested in millimetric laser ranging 
experiments in the Solar system. 



VII. CONCLUSION 

Laser ranging and time transfer in solar system could 
be able to reveal smoking gun arguments against GR by 
measuring accurately deviation from GR in well-suited 
experiments such as LATOR [16], TIPO [17], ASTROD 
[TS], ODYSSEY QI] or SAGAS [2D]. Such experiments 
will require distance measurements at millimetric level. 



As emphasized in pQ 



terms in gij then have to be 



taken into account. But only monopolar terms have to be 
considered in the c part of the metric. Thus, definition 
of multipolar moments is required at the c -2 level only. 
Such definitions are available for GR and TS theories 
[TP] [14] and are in progress in the PPN formalism [11] . 
Moreover, since the Solar system is strongly hierarchized, 
only the Solar monopolar term has to be considered in 
the c -4 part of the metric. 

Hence, from numerical considerations that spring from 
present constraints on the post-Newtonian parameters 
and from the Solar system specificities, the 2PN/RM 
metric (51)-(58) is sufficient for the next generation of 



experiments dealing with propagation of light in the So- 
lar system. 
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Appendix A: Effects on the metric of a non-rigid 
rotation 

Recent results coming from Solar seismology suggest 
that the Solar core rotates faster than the external lay- 
ers [37-29 . In the following toy model we consider the 
mass density as spherical (p(x) = p{r)). As usual in So- 
lar models |25j , we make the distinction between three 
main regions : the core region (r e [0, Rc ~ 0,2i?s], 
where Rc is the core radius), the radiative region (r G 
[Rc,Rr ~ 0,7i?s], where Rr is the tachoclyne radius) 
and the convection region (r E [Rr, Rs], where R$ is the 
Sun radius). f_ being the angular velocity, we modelize 
the differential rotation as follows 

Q(x) = fl(r)k z , (Al) 

with 

S)(r) = Si c (9) n (-£] (A2) 

where 

IL(x) = 1VsG[0;1] 

= everywhere else, 

and where we modelize the differential rotation by a sim- 
ple model in accordance with the usual model |26j 

n B (9) =n B (l + e B cos 2 8), (A3) 

where Q B and e B are constants, B being any of the three 
previous regions (C, R or D). 

Solar seismology suggests that the radiative region ro- 
tates as a solid - meaning £Ir(9) = Qr (ie. €r = 0). 
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However, since we are interested in testing solar seismol- 
ogy results, let us relax this assumption. Let us write 



W l s (t, x) = Wh{t, x) + W R (t, x) + W* D (t, x), (A4) 

4 p(R) 



with 

W B (t,x) =G [ d 3 X (n B (R) 



x X 



\2-X\ 



(A5) 



Wg being the Sun spin part of Wi in (52). Then the 
solution writes 



(q b x x 
Where 



(A6) 



1 , e B \M? , ^Mf (5co ^_ i: 



15 J H 



35 r 5 



M% = M N (0,R c ), M* = M n (Rc,Rr), 

Mjy = M N (R R ,R 8 ), 



with 



M N (X,Y)= [ R N+2 p{R)dR. (A7) 
Jx 

First, note that a faster rotation of a rigid core will 
modify the value of the total angular momentum only. 



However, this value could be affected by a differential 
rotation as well. 



But, one also may have to consider a term like the last 
term of the r.h.s. 



of ( A6 ) in ( 56 ) , in order to measure 
possible weak effects due to differential rotations of the 
different stages of the Sun - then, giving a characteris- 
tic way to put constraints on such differential rotations, 
which will be independent of Solar seismology and neu- 
trino detection results. In what follows, we will refer to 
this term as the post-Lense-Thirring term. 



As an illustration, let us consider the following toy 
model. Assume (1) the density decreases linearly with 
the distance to the center of the Sun (2) the differential 
rotation is independent of the distance from the center. 
Now consider that the photons - used for the time trans- 
fer pQ - graze the Sun (ie. b = aRs, where b is the 
impact parameter and a{> 1) a parameter ideally close 
to 1). Then the effect of the post-Lense-Thirring term is 
about a~ 2 es/ll times the usual Lense-Thirring effect. 



However, realistic models of the Sun and its rotation 
are expected to substantially decrease the value of this 
possible effect. But, this point should be verified and a 
specific study, that consider different realistic models of 
the Sun, should be done in order to clarify this issue. 



